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ABSTRACT 

A convergent  aeries  representation  for  the  internal  partition 
function  of  polyelectronic  atoms  is  obtained  by  assuming  a covolume 
equation  of  state  for  the  gas  as  previously  applied  by  Fermi  and  Urey 
to  the  hydrogen  atom. 

The  present  investigation  is  limited  to  those  cases  wherein 
only  extr&nuclear  electronic  excitation  occurs.  The  contribution  of 
these  electronic  states  to  the  thermodynamic  functions  is  obtained 
from  an  acceptable  approximation  to  the  sum  of  the  convergent  sex  - 
ies  for  the  partition  function. 

It  is  shown  that  at  relatively  low  temperatures  (3000°K),  the 
customary  method  of  evaluating  the  Internal  partition  function  (based 
on  the  assumption  of  an  ideal  gas)  agrees  to  within  a few  percent  with 
the  results  obtained  from  the  covolume  treatment.  However,  at  high* 
er  temperatures  the  increase  in  size  of  the  excited  atoms,  along  with 
the  appearance  of  charged  particles  produced  by  ionization,  render 
the  ideal  gas  treatment  inadequate.  Since  the  interaction  potentials 
of  charged  particles  are  not  known  in  general,  an  approximate  pro* 
cedure,  which  neglects  these  interactions,  is  suggested  for  analyzing 
a system  wherein  ions  and  free  electrons  constitute  a small  fraction 
of  the  total  population.  This  procedure  should  be  useful  for  treating 
gaseous  mixtures  to  temperatures  of  about  10,000°K. 
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I.  INTRODUCTION  AND  SUMMARY 


If  the  temperature  of  an  atomic  gas  is  sufficiently  high  to  pro- 
duce appreciable  equilibrium  concentrations  of  atoms  in  excited  energy 
states,  then  an  accurate  calculation  of  the  thermodynamic  functions  of 
the  system,  including  contributions  from  these  higher  energy  states, 
becomes  difficult  because  of  the  apparent  divergence  of  the  partition 
function.  Practical  examples  of  high  temperature  systems  of  this  type 
arc  encountered  in  researches  on  stellar  atmospheres  and  in  shock 
waves.  This  problem  is  also  of  interest  in  connection  with  the  appli- 
cation of  high-performance  propellants  and  nuclear  energy  sources  to 
jet  propulsion  devices. 

It  has  been  known  for  some  time  that  the  internal  partition 
functions  oi  ideal  monatomic  gases  diverge  if  all  of  the  possible  energy 
states  are  included.  The  terms  of  the  infinite  series  which  define  the 
partition  function  consist  of  the  product  of  the  statistical  weight  and  of 
an  exponential  factor  involving  the  negative  of  the  excitation  energy  of 
each  quantum  state  divided  by  k.T  . This  exponential  factor  for  the 
higher  states  approaches  a finite  limit,  whereas  the  statistical  weight 
increases  indefinitely.  Thus  although  the  first  few  terms  of  the  par- 
tition function  may  appear  to  be  converging,  the  higher  terms  gradually 
Increase  and  become  arbitrarily  large.  The  series  therefore  diverges. 
However,  it  is  known  that  at  relatively  low  temperatures,  the  first  few 
terms  of  this  series  yield  a good  approximation  to  the  internal  parti- 
tion function. 

The  principal  resuits  obtained  in  the  present  analysis  are: 

(a)  a convergent  representation  is  obtained  for  the  partition  function 
of  polyelectronic  atoms;  (b)  relations  for  the  thermodynamic  functions 
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of  un-ionised  gases  at  very  high  temperatures  are  derived  from  this 
expression;  and  (c)  the  limit  of  validity  of  the  customary  procedure 
is  determined  by  comparing  results  for  the  thermodynamic  functions 
based  on  that  method  with  the  results  obtained  from  the  convergent 
series.  The  convergent  series  representation  of  the  internal  parti- 
tion function  for  all  atoms  is  obtained  by  selecting  a covolume  equa- 
tion of  state  for  the  gaseous  mixture.  This  procedure  was  used  prev- 
iously by  Fermi^)  and  Urey^  in  their  investigations  on  the  hydrogen 
atom.  Computations  based  on  the  present  method  show  that  the  co- 
volume correction  is  the  major  deviation  from  ideal  gas  behavior 
for  gaseous  mixtures  at  relatively  low  temperatures  (3000°K)  at  mod- 
erate and  at  elevated  pressures.  However,  it  is  found  that  this  cor- 
rection is  small,  and  that  the  convergent  solution  to  the  partition 
function  yields  results  which  agree  to  within  a few  percent  with  the 
results  obtained  for  an  ideal  gas.  Thus  it  would  appear  that  the  cus- 
tomary procedure  for  computing  internal  partition  functions  at  low 
temperatures  is  entirely  adequate  for  engineering  calculations.  At 
higher  temperatures  ions  and  free  electrons  appear  in  appreciable 
concentrations,  and  proper  account  must  be  taken  of  the  increase  in 
the  number  of  particles  although  the  interaction  of  charged  particles 
in  the  mixture  can  be  neglected  in  first  approximation* . By  applying 
the  present  method,  one  can  determine  the  relative  importance  of  the 
covolume  and  ionization  corrections.  Sample  calculatiohs  for  hydro- 
gen, lithium  and  nitrogen  show  that  at  temperatures  above  about 
7000°K  ionization  can  rot  be  neglected  even  at  pressures  as  high  as 
20  atmos. 

* This  effect  may  be  accounted  for  by  using,  for  example,  the  Debye- 

Htickel  theory  as  applied  by  Williamson  in  the  case  ol  hydrogen. 

Cf.  Astrophysical  J.  103  139  (1946). 


The  procedure  developed  for  the  calculation  of  the rrr.odyna.mlc 
functions  of  neutral  atoms  remains  valid,  in  principle,  even  under 
conditions  in  which  appreciable  ionization  occurs.  However,  before 
applications  can  be  made  to  real  systems,  it  will  be  necessary  to  in- 
vestigate in  some  detail  (a)  small  perturbations  of  the  calculations 
for  neutral  atoms  and  (b)  interaction  terms  corresponding  to  binary 
collisions  between  neutral  atoms  and  charged  particles.  But,  these 
calculations  are  extremely  difficult  to  carry  out,  because  the  inter- 
action potentials  of  charged  particles  are  generally  not  known. 

It  is  expected  that  these  interaction  terms  will  yield  small 
corrections  on  the  thermodynamic  functions  of  systems  wherein  ion- 
ization is  of  the  order  of  10  percent  (temperatures  at  about  10,  000°K).  ' 
Consequently,  a first  approximation  to  these  functions  may  be  obtained 
by  treating  the  system  with  the  methods  of  classical  solution  thermo- 
dynamics and  neglecting  the  interaction  between  the  neutral  atoms 
and  the  charged  particles.  The  present  method  may  be  utilized  to 
compute  the  thermodynamic  functions  of  the  neutral  atoms,  and  the 
contributions  of  the  ions  and  free  electrons  may  be  calculated  by  ua» 
ing,  for  example,  the  Debye -Hiickel  theory. 

The  problem  of  evaluating  partition  functions  may  be  attacked 
by  using  either  the  methods  of  statistical  mechanics  or  the  methods  of 
thermodynamics.  The  former  approach  was  applied  by  R.  H.  Fowler^ 
and  by  Planck*  to  the  investigation  of  stellar  atmospheres.  The  latter 
approach  was  applied  by  Fermi^  and  by  Urey^.  The  advantages  of 
the  statistical  methods  can  bo  utilized  only  if  a realistic  description 
of  the  interaction  potentials  for  atoms  in  their  ground  states  and  in 
*C i.  Reference  3,  p.  353. 
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excitei  electronic  state3  can  be  given.  J:ince  these  interaction  poten- 
tials are  generally  not  .no'.vn,  the  statistical  methods  offer  no  obvious 
advantages  over  the  simpler  thermodynamic  methods,  which  involve 
the  implicit  assumptions  that  the  atoms  are  rigid  elastic  spheres  in- 
teracting according  to  a square-well  potential  profile. 

The  present  study  utilizes  the  methods  of  statistical  thermo- 
dynamics for  the  evaluation  of  internal  partition  functions  for  atomic 
gases  distributed  at  equilibrium  in  the  accessible  electronic  energy 
states.  It  is  assumed  that  the:  presence  of  atoms  in  excited  energy 
states  leads  to  non-ideal  g&s  behavior  only  insofar  as  a covolume  cor- 
rection is  required  to  the  equation  of  state.  As  is  to  be  expected  from 
the  use  of  a covolume  correction,  the  necessary  weighting  function  to 
provide  convergence  of  the  series  defining  the  internal  partition  func- 
tion is  dependent  upon  the  effective  atomic  volumes  of  the  excited 
electronic  states.  This  conclusion  was  also  stated  by  Urey  and  by 
Fermi.  It  is  shown  in  the  present  studies  that  the  exact  expression 
for  the  weighting  f unction  can  be  obtained  along  with  the  solution  for 
the  total  partition  function.  The  reduction  of  the  present  solution  to 
the  classical  value  of  the  partition  function  for  an  ideal  gas  is  readily 

demonstrated. 

Assuming  the  validity  of  a covolume  equation  of  state,  it  is 
clear  that  the  calculation  of  the  thermodynamic  properties  of  gas- 
eous systems  depends  upon  the  evaluation  of  the  weighting  function. 

It  is  found  that  the  weighting  function  is  determined  by  the  solution 
of  three  simultaneous  transcendental  equations  wherein  the  unknowns 
are  dependent  upon  the  quantum  energy  levels  and  upon  the  thermo- 
dynamic parameters  of  state,  viz,  , the  pressure  and  temperature. 
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A solution  to  the  set  of  rirrultaneous  equations  is  obtained  by  intro- 
ducin'; a closed  form  expression  for  the  internal  partition  function. 

The  dependence  of  the  weighting  function  on  the  atomic  ’ml- 
■irr.es  has  been  Indicated.  T.vo  methods  for  estimating  these  volumes 
are  investigated.  The  first,  an  approximate  techniciue,  is  the  method 

(4,0 

o£  screening  constants;  the  second  uses  empirical  wave  functions  ob- 
tained by  the  method  of  Hartree^  from  known  energy  levels.  The 
latter,  although  more  reliable,  is  extremely  laborious.  It  is  shown 
by  means  of  actual  calculations  on  the  potassium  atom,  that  the 
screening  constant  treatment  yields  atomic  volumes  which  are  in 
good  agreement  with  those  obtained  by  the  Hartree  method,  especially 
for  the  higher  energy  states.  For  these  states,  the  effective  core 
charge  reduces  to  one  so  that  an  atom  acquires  essentially  a hydrogen- 
like structure.  It  is  believed  that  the  application  of  the  screening  con- 
stant method  to  the  lower  energy  states,  and  the  use  of  the  hydrogen- 
like approximation  for  the  higher  states,  provides  a useful  teclmique 
for  evaluating  atomic  volumes. 

Application  of  the  proposed  method  to  the  determination  of  the 
thermodynamic  functions  of  pure  atomic  gases  at  high  temperature  is 
illustrated  by  Gample  calculations  on  the  un-ionized  lithium  atom. 

The  calculations  have  been  performed  for  pressures  from  1 atm.  to 
100  atm.  , and  for  temperatures  between  5000°K  and  20,  003°K.  The 
results  obtained  for  100  atm.  pressure  are  compared  with  the  classical 
solution,  and  it  is  found  that  in  the  temperature  range  wherein  ioniza- 
tion is  small  the  ideal  gas  results  differ  by  a few  percent  from  the  re- 
sults based  on  a covolume  equation  of  state. 


H.  THERMODYNAMIC  DERIVATION  OF  A GENERAL 


. EXPRESSION  FOR  THE  PARTITION  FUNCTION 

A.  Co  volume  Equation  of  State 

At  high  temperatures  appreciable  numbers  of  atoms  will 
occupy  higher  energy  states  for  which  effective  atomic  volumes 
are  n . longer  negligibly  umail*  hence,  the  need  for  a covolume 
correction  to  the  equation  %J  *i*.t e.  It  is  apparent  that  the  atomic 
volumes  increase  with  increasing  energy  levels.  Thus,  the  avail* 
able  free  space  in  the  gas  is  decreased  thereby  increasing  the 
probability  for  collision.  Several  investigators  have  examined 
the  problem  of  collisions  between  rigid  elastic  spheres,  and  it  is 
agreed  that  on  the  basis  of  such  a model,  the  covolume  correction 
b is  proportional  to  the  average  effective  volume  bm  of  a binary 
collision  (the  gas  densities  of  interest  here  are  too  low  to  yield 
appreciable  probabilities  for  collisions  involving  more  than  two 
atoms).  First  order  calculations  for  the  proportionality  factor 
indicate  |r  to  be  the  correct  value^.  In  this  case  N is  the  total 
number  of  atoms  present  in  the  mixture  of  total  volume  V , pres- 
sure P and  temperature  T . The  corresponding  equation  of  state  is 

P(v-b)  = P(V-  = NlcT  (1) 

Following  Ferml^  and  Urey^'  we  consider  a Bystem  which 
coneists  of  only  one  chemical  component  and  regard  the  various  ex- 
cited electronic  states  of  the  atom  as  separate  species.  The  re- 
sulting system  is  composed,  therefore,  of  a mixture  of  gases  of  the 
same  atomic  weight.  Of  the  N atoms  present  in  this  mixture,  Nj 
will  occupy  the  internal  energy  state  0JL.  Evidently 
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7]Nl  = N (2) 

T 

The  average  effective  collision  volume  bm  may  be  evaluated 
in  terms  of  the  individual  atomic  volumes.  Consider  therefore  the 
collision  of  an  atom  of  species  i with  one  of  species  j.  If  we  denote 
by  and  the  radii  of  these  atoms , then  the  effective  collision 
radius  ffj:  for  rigid  elastic  spheres  is 


^ - n.  + ^ 


where  b • b * volume  of  atom  of  type  L » i • The  effective  collision 
1 J o 

volume  bu  i*  equal  to  that  of  a sphere  of  radius  f-  . Thus 

0 .i  S 


I 1 3 

* L ^ . -3 


br4f(V5-  (bzbz 


An  average  effective  collision  volume  b for  the  system  may  be  de- 

m 

fined  by 


b_s  Z 


[total  number  of  colllsions)x(corresponding  collision  volume) 

(total  number  of  collisions) 


V 

The  total  number  of  collisions  between  atoms  of  types  t and  j is  given 

by  the  number  of  combinations  of  N|;  and  N things  taken  two  at  a time 

L ( 

such  that  one  of  each  is  present  in  every  set.  Thus  the  total  number 
of  collisions  is 

T.  N;  N; 


bm  = --ty 

ZNtNj 


7 


I 


But  from  Equation  (2),^NL-N  ; therefore,  NlN;=Ntand 

li 

* b"=  N2r.NiNi(bi!+bj')5 


!3) 


Hence  the  covolume  correction  b In  terms  of  the  atomic  volume  is 

b~2N]CNiNj  (bi  + b^>)  (4) 

This  result  was  indicated  by  Fermi,  ^ but  without  proof. 

The  generalized  expression  for  the  covolume  correction  b 
based  on  a two-body  collision  mechanism  is  reexamined  in  Sec.  VI  C in 
greater  detail.  A second  relation,  more  suitable  for  equation  of 
state  calculations,  is  readily  obtained  from  this  form. 


B.  The  Total  Partition  Function 

It  is  the  purpose  of  the  present  section  to  derive  an  expres- 
sion for  the  total  partition  function  of  a gaseous  system  of  chemi- 
cally identical  atoms  occupying  various  electronic  energy  states 
by  using  the  methods  of  classical  solution  thermodynamics 
(cf.  Appendix  A).  Statistical  relations  will  be  used  only  to  assist 
in  identifying  and  simplifying  the  results.  The  final  expression  for 
the  partition  function  of  imperfect  gases  is  shown  to  be  of  familiar 
form,  differing  from  the  classical  solutions  only  through  an  addi- 
tional factor  which  represents  a covolume  correction.  Furthermore,, 
this  correction  factor  is,  In  zeroth  approximation,  identical  with 
the  result  of  Fermi,  and  is  also  of  the  same  form  as  the  CO  r recti  on 
factor  obtained  by  Urey  for  a mixture  of  hydrogen  atoms  obeying  the 
covolume  equation  of  state. 
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An  expression  for  the  total  partition  function  can  be  derived 
from  the  definition  of  the  chemical  potential.  .Let  us  consider  the 
chemical  potential yU-  for  an  atom  of  species  i in  a gaseous  mix- 
ture composed  of  atoms  which  are  in  different  energy  levels.  Then 
(Cf.  Appendix  A), 

IcTinMi  (A-3) 

where  u°  is  the  chemical  potential  of  l in  the  puru  stats,  )o  rep- 
resents the  Boltzmann  constant  and  is  evidently  the  mole  frac- 

N 

O Q 

tion  of  i . If  h:  anr*  A,  denote  the  enthalpy  and  entropy  per  atom  of 
pure  component  i at  P and  T in  a volume  VL°  occupied  by  NL  atoms, 
then  we  may  write 


o 


(5) 


Note  that  the  superscript  o denotes  a pure  component. 

The  total  energy  per  atom  may  be  separated,  in  good  ap- 
proximation, into  translational  and  internal  components.  Thus 


t[  = | k.T  + 


(6) 


The  enthalpy  for  NL  atoms  of  species  l in  the  pure  state  at  pressure 

o 

P , temperature  T , and  volume  is 

+ PVL° 

Hence,  the  enthalpy  per  atom  is 


_ Hf_  c*  , 

- N ~ + 


(7) 
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vvhere  the  total  internal  energy  E is 

e - T,  e ‘ = X Nte;  (A-21-, 

*•  l 

and  the  total  volume  is  given  by 

V * X!  vf  (A —17) 

L 

Substitution  of  Equations  (21.  (4)  and  (A-17)  into  (1)  yields 


y I PV*-  5i!iT'N1(b;5+b;>)3_  N.  IcT 

AJ  l 2 N ■ 4 ' L J ' L 

1 L 


t = 0 


This  equation  must  hold  for  a syste  IT1  C C Tltcl  ining  any  arbitrary  rai.r.. 
ber  of  atoms.  The  only  physically  acceptable  solution  is 

VM.f  h1 

Nl  2N 


-VL#a=£2Ni(bi’  + kj3)3+  kT 


(5) 


I 

It  is  evident  from  Equation  (4)  that 

J-EN,(bhtiV-  SKI 

2N  Y <T  L J ' dNL 
If  E quations  (o)  and  (9)  are  used  in  Equation  {7)  along  with,  if quatior 
(6),  then 


(9) 


h 


° = Sl  j i u).,  , p <|b 


Hence,  from  Equation  (5), 


^-IlcT+^  + Pjb  -U* 


and  from  Equation  { A - 3 ) 

In  5 i Mi  _ kti  _ P 

N 2 IcT  IcT  IcT^N^lc 

At  equilibrium «s  constant.  Therefore,  solving  for  N-  and 
summing  over  all  states,  -,ve  obtain 
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£ ^ = N - (N  e ?)  //U  y Ai K ATe  " £ ^ 

I ’ i. 

If  v/e  define  the  total  partition  function  Q In  the  visual  way,  ^ 


{1C) 


then 


^ / -fv  V A'A.  -"i/kT-  SL 

Q = (Ne  )2]e  e 

l 

-0. 

where  the  weighting  function  e <■  has  been  introduced,  and 


{ID 


a = P_  M) 

yi~  u aNL 


(12) 


3N; 


jj  or  ideal  gases  0:  is  evidently  equal  to  zero.  Since  the  factor 

Q 

increases  rapidly  with  increasing  energy  level,  e 1 will  cause  a 
rapid  decrease  of  the  numerical  values  of  the  terms  of  the  su.nrna- 
tion  as  1/Jj  is  increased. 

The  entropy  per  atom  -4°  may  be  separated,  in  good  approxi- 
mation, into  a translational  and  an  internal  component.  Thus 


'int 


n 


The  translational  component  (<1°)  is  independent  of  the  internal 


't^on-S 


(12) 


energy  states  It  i3  given  by  the  Sac  kur- Tetrode  relation. 

Tor  an  imperfect  gas,  the  appropriate  expression  for  the  transla- 
tional ent  ropy  per  atom  is 


(A,°  ) =A  In 

v L trails  trans 


jdJrmicTfe^V-b)1 

Nh3 


(14) 
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where  m is  the  mass  per  particle  and  h represents  Planck's  con- 
stant. The  contribution  to  4°  from  the  internal  energy  states  (4°) 
may  be  evaluated  by  means  of  the  following  argument:  It  is  known 
from  statistical  mechanics  that  LiZ. ) may  be  defined  by  the  rela- 

1 Lnc 


(A?). 

l int 


(15) 


where  XiX  denotes  the  number  of  quantum  states  available  to  the  sys- 
tem. In  the  case  of  a system  consisting  of  the  pure  component  l 
which  has  available  to  it  only  one  quantum  state  Itk  , -UP  is  evidently 
the  statistical  weight  (j.  (degeneracy)  of  that  quantum  state;  hence, 


(16) 


and  it  follows  from  Equations  (13)  to  (16)  that 


(2 jt  m IcT)2  ez(v-b)gi 

Nh3 


(17) 


Substitution  of  Equation  (17)  into  Equation  (11)  yields 

q _ (gjrmk.T)f(V-b)yq  e"  % AT  ' (io> 

h3  4^3. 

where  the  factor  before  the  summation  sign  evidently  represents  the 
translational  partition  function  for  a gas  obeying  the  covolume  equa- 
tion of  state.  One  can  readily  identify  the  complete  expression  (13) 
with  the  familiar  approximation 


Q = Q,  'Q.  , 

trans  Lnt 


(W) 


* Cf.  Ref.  10,  p.  92. 
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whers 


and 


q (2  jt  mlcl" ) 2 (V-  b) 

^rar.s  Pj 


-<x)-./k.T  - 9. 


(20) 


(21) 


The  expression  for  the  internal  partition  function  Q.  . ia  of  the  usual 

int 


-0. 


Q = (2Jimk.T)lv  Vg  e~ 


form  except  for  the  factor  e 1 which  appears  as  a direct  consequence 
of  the  covolume  correction  introduced  into  the  equation  of  state.  It  is 

of  interest  to  note  that  as  -+-0>  Equation  (21)  reduces  to  the  total 

ai\iL 

partition  function  of  a perfect  gas: 

■ - \ / - as-.  /Alt 

H3  1 

It  is  well  laiown  that  this  expression  diverges  if  enough  terms 
are  considered  in  the  summation.  The  customary  practice  has  been 
to  compute  Q.  , by  using  only  the  first  few  terms  of  this  relation  and 

LMt 

to  assume  that  the  higher  energy  states  make  negligible  contributions. 
At  low  temperatures  this  procedure  yields  acceptable  results.  For 
these  cases  the  comparison  of  results  obtained  from  such  an 
approximation  with  the  correct  values  obtained  from  Equation  (21) 
shown  agreement  to  within  one  percent.  It  is  of  interest  to  note  that 
for  temperatures  in  the  order  of  5000°K,  the  first  term  alone  from 
either  of  the  above  expressions  for  Q yields  results  which  are 
accurate  to  within  a few  percent. 

The  functional  form  of  the  covol  dine  correction  given  in  Equation 
(12)  is  identical  with  the  relation  obtained  previously  by  Formi^^  as 


the  result  of  a less  complete  analysis  than  haa  been  given  here.  It  is 
evident  from  Equation  (20)  that  a covolume  correction  occurs 


also  In  the  expression  for  the  translational  partition  function.  How- 
ever, this  term  is  oi  no  consequence  as  regards  the  calculation  of 
internal  partition  functions.  Equation  (21)  will  be  referred  to  here- 
after ac  Fermi's  equation  for  the  internal  partition  function  of  a gas 
mixture  obeying  the  covolume  equation  of  state  or,  more  briefly, 
as  Fermi's  equation.  It  will  serve  as  the  basic  relation  for  the  de- 
termination of  thermodynamic  functions  of  atomic  gases  at  elevated 
temperatures. 
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ni.  QUANTUM  STATES  FOR  POEYEEECTRONIC  ATOMS 
A . Excited  Electronic  States 

The  inten  al  partition  function  as  given  by  the  Fermi  Equation 
is  defined  as  a infinite  series  which  extends  over  all  internal  states 
i,  beginning  with  the  ground  state  1=1.  Each  of  these  states  i de- 
notes a distinct  quantum  state  of  the  atom.  In  the  present  analysis 
we  shall  confine  our  attention  to  a consideration  of  the  states  of  the 
extra-nuclear  electrons,  and  shall  neglect  the  quantum  states  of  the 
nucleus.  The  rotational  and  vibrational  excitation  energies  of  the 

nucleus,  even  for  the  heaviest  elements,  ^ ^ are  known  to  be  at  least 

4 

of  the  order  of  10  e.  v.  These  energies  are  several  orders  of  mag- 
nitude larger  than  the  highest  degree  of  extranuclear  excitation  which 
can  be  expected  at  the  temperatures  of  interest.  At  most  we  shall  be 
concerned  with  the  first  few  ionization  stages  of  the  atom.  Consequent- 
ly, the  index  l is  used  to  represent  the  excited  electronic  3tates  of 
the  atom  and,  in  fact,  may  be  replaced  by  an  appropriate  set  of  quan- 
tum numbers.  Thus 

i = L (n,  L , 5 ) 

where  n is  the  principal  quantum  number,  L the  quantum  number 
measuring  the  total  angular  momentum  vector  for  the  atom,  and  3 , 
the  spin  quantum  number  for  the  atom. 

The  excitation  energies  ur.  are  given  in  the  usual  way,  viz.  , 

to-.  = W,  - W- 

L 1 L 

Here  VV.  denotes  the  term  value  for  state  ( and  Wt  refers  to  the 
ground  state.  These  data  are  readily  available  from  spectroscopic 
measurements  (e.g,  Eandolt  - Born3tein  "Tabelisn1  ).  An  accurate 
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cclculation  of  Q.  , requires  c.  complete  compilation  of  thee*-  term 
Int 

values.  Although  there  exist  infinitely  many  term  values,  the  usual 
practice  in  tabulating  these  data  is  to  list  the  values  Wt  up  to  perhaps 
n = 10.  For  temperatures  below  3000°K,  it  has  been  customary  to 
compute  Q.  by  utilizing  only  the  first  few  terms  in  liquation  (21). 

At  higher  temperatures,  it  is  neceaoary  to  take  more  terms. 

B.  Hydrogen- L.ike  Approximation  for  Higher  Electronic  States 

It  has  been  noted  that  the  higher  electronic  term  values  may  not 
have  been  determined  experimentally;  nevertheless , it  is  possible  to 
predict  these  term  values  by  assuming  a hydrogen-like  atomic 
configuration.  Thi3  approximation  can  be  justified  on  the  basis  of 
simple  physical  arguments.  As  the  excited  electron  occupies  higher 
energy  states,  the  net  effect  of  the  core  structure  (the  nucleus  and 
remaining  electrons)  on  thi3  electron  is  much  like  the  effect  of  the 
hydrogen  nucleus  on  its  electron.  The  inner  electrons  shield  the  excited 
electron  from  the  full  strength  of  the  nuclear  field.  In  fact,  as  the  atom 
approaches  the  ionized  state,  the  effective  core  charge  approaches  1, 
and  the  structure  becomes  hydrogen-like.  In  practice  this  limiting 
condition  is  assumed  to  exist  much  earlier.  Consequently,  one  can. 
estimate  approximately  the  remaining  term  values  by  using  the  well- 
known  Bohr  hydrogen  model. 

It  appears  appropriate,  therefore,  to  separate  the  internal 
partition  function  into  two  parts.  The  first  part,  a finite  sum,  accounts 
for  the  contribution  of  the  known,  spectroscopically  determined,  tern- 
values.  For  lev/  temperatures  this  part  makes  the  major  contribution 
to  Q . The  remaining  portion  of  the  internal  partition  function  ia 
determined  by  a hydrogen-like  approximation,  and  ia  represented  by 


an  infinite  series;  thus 
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t 


> 


-w,  /kj 


v w/iCT-el  ’ 

Z-9ie  +Z3,e 


W./JcT-0, 


(22) 


T ~ l 

...  3ix' 

L— 1 i=^+l 

The  separation  is  made  at  some  value  of  l—/  3uch  that  W is  de- 

L 

pendent,  in  first  approximation,  only  on  the  principal  quantum  num- 
ber n , that  istosay.it  is  assumed  that  for  given  n , the  term  values 
differ  by  negligible  amounts  for  the  allowed  values  of  L and  S . 

This  assumption  is  well  substantiated  by  the  spectroscopic  data.  An 
inspection  of  these  data  shows  that  for  a particular  atom,  there 
exists  some  value  of  n=a  , beyond  which  the  L and  S variations  in 
the  term  values  may  be  neglected.  It  follows  therefore  that  for 


W[  ~ wn 

In  accordance  with  the  hydrogen-like  approximation.,  the  index  £ will 
be  replaced  by  the  pair  n and  £ , the  principal  and  azimuthal  quantum 
numbers  respectively.  Thus 

The  infinite  series  in  Equation  (22)  may  now  be  written  as 


n=a  1=0 


K,At  -9„t 


(23) 


The  summation  over  i is  replaced  by  the  double  sum  over  n and  £ . 
Although  the  description  of  the  energy  states  by  the  quantum  numbers 
( n,  l ) is  a customary  one,  the  above  formulation  is  primarily  a for- 
malism. Note  that  for  given  n,  the  summation  over  the  allowed  values 
of  £ 1b  required  first;  the  summation  over  n i3  carried  out  last. 
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I 


I 


•sr 


This  procedure  has  no  particular  physical  significance  and  has  been 
adopted  only  because  it  i3  conve  i no  correct  order  of  the 

absolute  magnitudes  of  the  energy  states  ie  known  to  be* 

14,  24, 2 p,  34,3p,  44, 3d,  4 p,  54. 4 d,  5p,  4f,  <54  ... 

The  procedure  indicated  in  liquation  (23)  may  be  used,  because 
the  series  is  absolutely  convergent. 

If  the  hydrogen-like  approximation  is  applied  to  the  calculation 
of  0£n  and  ^ , then  one  can  readily  compute  appropriate 
expressions  for  each.  The  statistical  weight  0^  for  a hydrogen-like 
configuration  is  given  by  the  product  of  the  statistical  weight  of  the 
core  and  the  weights  for  the  single  excited  electron.  If  we  denote 
the  core  degeneracy  by  0 , and  note  that  the  degeneracy  of  a single 
electron  for  all  values  of  l for  given  n is  ZIZl  + 1),  then  we  find 
that  the  statistical  weight  for  a given  state  Uh  is 


9ln~  29c(2^+  0 i24) 


'"he  term  values  Wj  for  a hydrogen-like  configuration  are  given  by 

W,  ^ Rohe  125) 

n2 

where  Ro  is  the  Rydberg  constant  (13.53  e.  v.  ),  and  c is  the 
velocity  of  light.  The  substitution  of  Equations  (24)  and  (25)  into  the 
series  (23)  yields 


* Cf.  Pauling,  L. , "The  Nature  of  the  Chemical  Bond",  Cornell 
Univer  *ty  Press,  Ithaca,  2nd  '.edition  1940. 
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t 


The  state  parameter  od  is  equal  to  R0hc/k.T . 

The  internal  partition  function  may  be  given  now  in  terms  of 
n and  f. : 


Q = e 

hnt 


rw'ATi 


■A  » , p n-i  A 

Iq  + 29Ze^lW*-‘ 

1=1  *q  *“x  n=Q  l-o 


(26) 


where 


w.  At- 0 

= a i'  i 


% s 9te 


(27) 


C.  Ionization  and  .V.ultiplc  Joxcitation  of  Polyelectronic  Atoms 

The  present  study  is  limited  to  a consideration  of  gaseous 
systems  at  temperatures  and  pressures  for  which  the  chemical  species 
present  will  occupy,  in  appreciable  numbers,  only  the  electronic  states 
of  the  neutral  atom.  A system  of  this  kind  ie  encountered  in  several 
engineering  problems  which  are  of  current  interest,  viz.  , high-per- 
formance rocket  devices  and  very  strong  shock  waves.  Important  ex- 
ceptions occur  in  problems  of  stellar  atmospheres  where  extensive 
ionization  i3  present.  In  the  outer  regions  of  the  stars  (chromosphere) 
the  pressures  are  of  the  order  of  10  ^ atmospheres,  and  under  these 
conditions  the  constituent  elements  of  the  atmosphere  (e.  g.  Ca,  ila, 

Sr,  Si,  He.  . . . ) appear  in  the  fir3t  and  higher  stages  of  ionization,  i ‘al- 
though the  present  analysis  is  derived  for  systems  which  contain  only 
neutral  particles,  it  is  evident  that  similar  treatments  can  be  developed 
for  mixtures  of  atoms,  ions  and  free  electrons,  however,  the  results 
which  we  obtain  are  useful  for  the  description  ol  the  un-ionizea  atoms 
even  in  mixtures  which  contain  charged  particles.  If  the  ions  and  elec- 
trons constitute  a small  fraction  of  the  total  population,  then  one  can 


ft 
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ft 


obtain  an  approximate  treatment  by  neglecting  their  interaction 
with  the  neutral  atoms. 

The  applicability  of  the  present  method  to  a specific  problem 
depends  on  the  degree  of  ionization  which  occurs  in  the  system,  and 
the  proposed  2>rocedurc  is  valid  as  long  as  the  number  of  ionized 
atoms  represent  only  a few  percent  of  the  total  population.  The  de- 
gree of  ionization  for  a particular  component  in  a mixture  is  inde- 
pendent of  the  other  elements  present.  Thus  v/e  may  confine  our 
attention  to  a single  reaction,  which  may  be  represented  by  the  equa- 
tion 

X^X+-t-  e (28) 


X is  the  chemical  symbol  of  the  neutral  element  and  e represents 
the  electron  freed  in  the  ionization  process.  It  war.  pointed  out  by 

/ 14) 

Saha'  that  this  reaction  could  be  treated  by  standard  thermody- 
namic methods;  therefore,  the  concept  of  an  equilibrium  constant 
is  applicable.  L.et 

N0  = N - Nl  A N = 1 

(2  ) 

N,  = Ne  T]_=  N + Nl 


where  AN  denotes  the  net  increase  in  the  number  of  particles  produced 
by  ionization  as  shown  by  liquation  (28),  and  T|  the  total  number  of 
particles  present  at  equilibrium.  N represents  the  initial  number 
of  atoms  of  X , N0  the  number  of  neutral  atoms  at  equilibrium,  Nt 
the  number  of  ionized  atoms , and  Ne  the  number  of  electrons.  The 
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eq’j.i.lihriurri  constants  arc  related  by 


D 

KP=  ^N(f\) 

and  K,m=  NiNe//Na  . IJvidcntl; 


PN, 


-r  N -N, 

If  denote  Of.  = N. /TnI  * then 


= ( d/p 


l 


X = i-X 


1+<P/P, 


(3d) 


Therefore,  the  degree  of  ioniration  is 


N,  _ x 

N0  1-^ 


(•31) 


The  equilibriurv  constant  K in  obtained  readily  from  the 
familiar  expression 


which  corresponds  to  the  reaction 


Z>,A,^X>A 

1 k. 


and 


ibx-x>A'aAF”-AH‘-TAS' 

l k j 


1 


pi  O 

The  quantity  f-  . denotes  the  Gibbs  free  energy  of  the  pure  component 
A,  at  T and  P , and  is  related  to  the  free  energy  FI  A in  the  ;nb;= 


tore  by  an  expression  analogous  to  (G -3).  Note  that  this  formulation 
differs  from  the  customary  one  in  that  the  standard  properties  are  those 


9 
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of  the  pure  component  at  T and  P rather  than  the  properties  of 
A,  in  the  mixture  measured  at  standard  temperature  and  pressure. 
If  one  expands  the  expression  for  A F , noting  that  the  individual 
atoms  have  both  translational  and  electronic  degrees  of  freedom. 


then  one  obtains  the  v/ell-known  Saha  equation: 


.(15) 


I 


1C  - 2 (Q„t)i  / 2 JT me  k r 'I z kTc  lV‘ /kT 

(Q.AA  > 


(32) 


Here  me  denotes  the  mass  of  the  electron,  is  the  ionization  po- 

tential for  reaction  (28),  and  (Q^nt)j  is  the  internal  partition  function 
for  the  j -fold  ionized  atom.  The  substitution  of  this  expression  for 

into  Equations  (3'o)  and  (31)  yields  the  ratio  Hi  . Unsdld^  has 

No 

constructed  a nomogram,  v/hich  simplifies  numerical  computations 
appreciably. 

v.  hen  collecting  the  spectroscopic  data  for  the  evaluation  of 
the  finite  sum  in  Equation  (2o),  it  i3  important  to  note  that  the  tabu- 
lated data  will  necessarily  include  only  observed  transitions.  In  the 
case  of  an  equilibrium  process  (which  is  of  interest  in  the  present 
analysis)  all  energy  states  are  available  to  the  system.  Accordingly, 
one  must  calculate  term  values  which  are  not  observed  spectroscop- 
ically. Furthermore,  multiple  excitations  (the  simultaneous  excita- 
tion of  two  or  more  electrons)  may  occur.  An  accurate  computation 
must,  of  course,  make  allowance  for  ail  of  these  possibilities.  Al- 
though the  available  stales  can  be  predicted  on  the  oasis  of  the  possible 


mm 

electronic  configurations  by  well-known  procedures,  ' the  corres- 


*C£.  Reference  15,  pp.  66-74. 


« 
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ponding  t erm  values  arc  somewhat  more  difficult  to  specify.  The 

(17) 

method  of  isoelcctronic  sequences  may  be  useful  in  this  connec- 
tion. The  term  values  of  unobserved  states  for  a given  atom  may 
be  predicted  by  examining  the  corresponding  states  of  other  atoms 
which  have  the  same  electronic  configuration.  Because  of  similar 
electronic  configurations,  the  same  states  are  assumed  to  be  ac- 
cessible to  each  atom.  The  only  difference  between  these  various 
atoms  will  be  the  magnitudes  of  the  transition  energies  (term  values) 
because  of  the  different  nuclear  charges.  In  the  present  analysis  it 
was  found  convenient  to  plot  the  transition  energies  of  the  state  in 
question  as  a function  >£  the  effective  atomic  number  %eff  ' I^u8“ 

trative  plots  for  the  nitrogen  I configuration,  states  Si(l<f«24  2p?,  3d) 
and  D ( 14,24  2py,  are  presented  in  figure  1.  Evidently  one  or  two 
points  arc  adequate  for  establishing  z given  curve.  This  result  is 
produced  by  the  fact  that  the  transition  energy  is  proportional  to  the 
square  of  the  effective  field  strength. 
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IV.  CALCUtxATIOH  OF  THERMODTNAMIC  FUNCTIONS 


A.  Evaluation  of  W sighting  F unction 

The  derivation  of  the  weighting  function  prevented  in  Section 

-0 

II  showed  the  dependence  of  e 1 on  the  derivative  . The  defin- 
ition of  this  derivative  as  given  by  Equation  (9)  may  be  simplified  by 
introducing  the  Boltzmann  relation. 


N.  = e 

j 9i 


-wj/kT-ej+0! 


(33) 


In  addition  let  us  assume  that  the  atomic  volumes  b;  are  given  by 


1 


the  following  relation  ( /0 , aa  and  vj  are  defined  in  V A): 

br(^-)f 


(34) 


The  justification  for  this  assumption  will  be  presented  later.  The 
substitution  of  relations  (33)  and  (34)  into  (9)  yields 

3b  _ ao.'jrnjie^^  f q (V  +V  0/fcT-9i 

3Nt  ^ — 5 — /2g,N  tn'i  ^ 

But  it  is  known  from  statistical  mechanics*  that  the  internal  parti- 
tion function  may  be  defined  as 


' Q - &N 

N, 

and  from  the  Fermi  equation 


(35) 


-w,/kT  A w,/ict-9; 

Eae  1 '-e  M 

}-!  r‘  ‘ 


♦ Cf.  Ref,  12,  p.  116. 
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The  equation  for  may  be  written 


0, 


9i  ~ § (fsc)  ~ 


"I  ^ ( , , ?n  "i/v  ~wi 

A e '4»/'vj+v3je  J 


3 w;/jci-6. 


■i" 


oo 

XH, 

(-1  < 


j 


where  the  state  parameter  6 is  defined  by 


a=  (AxalXo\  P 
1 \ 3 k)C 


'ZVJ 

Equation  (36)  may  be  expanded 


(36) 


(37) 


0L-  Ae^1'  + ^J. 


2 '!’i  Ht  2 ^ v{ 

^n,+3vij=i_J+  .r*1 


P*  ’ &,  & 

* « d 


This  relation  cam  be  simplified  by  means  of  the  following  definitions: 


Thus 


=z-t% 

i-  1“ 

r-} 

III 

Lu 

r1  ' 

P*V’ +,“>;• 

+ H.  + T) 

(38) 


(39) 


And  note  that 


0i  = /6eVi5+/'<!+  W+'1) 


(40) 
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As  a convenience  in  terminology  we  shall  refer  to  yu.  , & and  Z as 
the  firsti  second  and  third  "moments"  respectively  of  the  internal 
partition  function. 

The  internal  partition  function  given  by  the  Fermi  equation 
may  be  written  in  terms  of  the  moments  yu  , £ and  T . 


Q « e W,/kJE 

j- 


<4i) 


Or,  in  the  notation  of  Equation  (26); 


e 2_^J)e 

1=0 


i 


It  is  assumed  here  that  Vn  is  a function  only  of  n » and  it  will  be  shown 
later  that  in  fact  V0=- n2  in  the  range  wherein  the  hydrogen-like  approxi- 
mation is  valid.  The  moments  yU  , & and  Z are  independent  of  l . Thus 
the  above  expression  reduces  to 


n-j 


where  the  well-known  relation  /^(2(  + i)=  n2  has  been  applied. 


l=o 

The  unknowns  Jj  , £ and  X defined  by  Equations  (38)  deter- 

-0 

mine  the  weighting  function  e 1 by  means  of  relation  (39).  This  equa- 
tion exhibits  the  direct  dependence  of  6j  on  the  atomic  volume  bj 
through  the  term  /3e9iVj3  (Cf.  Equation  (34)).  This  term  is 
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easentially  the  solution  obtained  by  both  Fermi  and  Urey  in  their 
investigations.  It  is  important  to  note  that  all  four  terms  in  dj  must 
be  included  in  the  evaluation  of  the  weighting  functi  on  In  fact,  for 
th  ' a e t^rrr.s  which  make  the  largest  contribution  to  Qi„i  > the 

effect  of  each  of  the  factors  involving  the  three  moments  are  equally 
important  and  of  the  same  order  of  magnitude  as  the  atomic  volume  term, 

. For  the  larger  values  >(  ft  ( T < 3000 °kl,  P>  1 atm, 

the  ma:or  contribution  to  the  internal  partition  function  is  from  the 
lower  energy  states,  as  is  well  known  The  relatively  email  magnitude 
of  e ^ cauces  the  rapid  convergence  of  the  series  (Cf.  Equation  (22)  ) 
so  that  the  higher  energy  states  make  in  all  negligibly  small  contri- 
butions to  The  function  0j  behaves  essentially  like  j6  , and 

its  value  for  the  lower  energy  states  is  small  relative  to  /icT 

(Cf  Equation  (59)  )•  Because  of  this  fact  its  exact  value  is  inconse- 

0i  1 

quential,  and  the  inclusion  of  only  the  volume  term  ^>e  \J,  yields  a 
suitable  order  of  magnitude  estimate.  Thus  the  approximation  suggested 
by  Uxey  of  using  only  this  term  in  0j  serves  to  demonstrate  the 
convergence  of  the  partition  function,  and  is  not  meant  to  yield  an 
accurate  computation  of  the  function.  For  the  larger  values  of  good 
estimates  of  can  be  obtained  readily  from  the  first  few  terms  of 

the  infinite  series  This  is  not  the  case  for  smaller  values  of  /S 
( T large,  p small).  Here  the  major  contribution  to  stems  from 
the  higher  energy  states  wherein  0^  is  of  the  order  of  the  energy 
term  Wj/k.T  . An  accurate  estimate  of  0^  is  therefore  required,  and 


one  cannot  omit  the  terms  involving  the  moments  jj 


• o 


and 


Obviously  the  Urey  approximation  is  no  longer  vaiid,  and  a more  exact 
solution  must  be  sought- 


~27 


It  is  evident  from  Equations  (38)  that  the  evaluation  of  the 
weighting  function  requires  the  solution  of  three  simultaneous 
transcendental  equations  for  the  quantities  yu  , b and  T.  It  is  noted 

Q 

that  each  of  these  relations  Involves  the  factor  y5e  1 through  the 

„ 0 

terms  . If  we  define  this  factor  p®  y3e  \ then 

i 


^U+yUv'p+M  + T) 


, — > j 2 

3ut  p K t vlyU  + Vi  $ + m'j « i ; therefore , 

A 

i - PU3+M2+^j  + t) 

When  P><  10 the  t term  in  the  bracket  above  is  predominant,  and 
the  entire  factor  is  well  approximated  by  px  • Thus 


1-yftT 


(43) 


At  lower  temperatures  (^>10*°),  the  terms  yuV£  and  ^ maybe  OC^) 
so  that  they  can  no  longer  be  neglected;  however,  in  these  cases  the 
entire  bracketed  term  is  so  small  that  it  may  be  entirely  neglected 
as  compared  to  1 . Thus  the  expression  (43)  is  applicable  to  the 
larger  values  of  ft  as  well  as  to  the  relatively  small  values  for  which 
it  was  specifically  derived. 

The  substitution  of  the  approximation  (43)  into  Equations  (40) 


and  (27)  yields 

3,e 


It  is  convenient  to  introduce  a quantity  q such  that 


(44) 


*28“ 


Then 


i' 

4 


i-Sv 


r1  J 


(45) 


Obviously 


so  that 


_ 3 (i-/3 1) 

/flcj 

- 2(L-£9M 

/6cj  ^ 

L/i-fi'i'f&i 


It  is  more  convenient  to  define 


(46) 


(4?) 


X(/J.4.t)-<ji  + a(l-^)|j-0  (48) 


A practical  method  for  solving  these  simultaneous)  equations  for  the 
unknowns  jj  , o And  T consists  of  a simple  iteration  procedure  based 
on  initial  estimates  for  these  quantities.  If  these  values  are  denoted 
by  \J0  , 80  and  ta  , then  the  functions  <J>  , X and  ^ may  be  expanded 


in  a Taylor  eerie*  about  the  point  (^0,&>,  0 . Thus,  for  example 


+ 4(|f )0A^+  (tt)0At  + •*  • -4<?* 


with  similar  expressions  for  the  other  two  functions.  If  only  the 
linear  terms  are  retained,  then  the  three  expansion*  reduce  to  three 
linear  simultaneous  equations  in  the  deviations  AyJ  , A & and  A't  . 
The  application  of  these  first  order  corrections  to  the  initial  esti- 
mate e yji0  , b0  and  T0  yield  a second,  more  accurate  solution.  The 
process  can  be  repeated  until  the  desired  degree  of  accuracy  is  ob- 
tained. 


B.  Internal  Partition  Function 


Evidently  a closed  form  expression  for  would  be  useful 
for  the  purpose  of  solving  Equations  (47).  Such  an  expression  for  C| 
may  be  obtained  by  deriving  a suitable  approximation  for  the  function 
, where 

X = 2Stn  e 


according  to  the  notation  of  Equation  (42).  Clearly,  this  function 

represents  the  individual  terms  of  Q.  for  the  higher  energy  states. 

mt 

If  one  plots  the  contributions  to  Q from  each  value  of  n (i.  e. 
the  terms  ^ ) u a function  of  n , the  area  under  the  function  will 
represent  Q . It  should  be  noted  that  the  lower  terms  in  the  finite 
sum  of  Equation  (42)  must  be  collected  for  given  n . Now  it  can  be 
shown  that  for  relatively  large  values  of  the  state  parameter  f})  , 
the  function  has  the  form  indicated  by  Diagram  A,  whereas  for 
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rela.ti.vely  small  value#  of  yS  , cj  behaves  essentially  as  shown  in 
Diagram  B. 


ft 

W 


The  step-like  nature  of  Q for  small  n represents  the  terms  of  the 
finite  sum  in  Equation  (42).  The  evaluation  of  the  partition  function 
for  the  cases  represented  by  Diagram  A la  a relatively  minor  prob- 
lem. Indeed,  it  is  of  the  type  which  one  encounters  at  the  lower 
temperatures.  For  these  casee  the  first  few  terms  , • • • • ^ 

yield  a sufficiently  accurate  estimate  of  G>.  , and  one  can  neglect 

int 

the  corrections  due  to  the  moments  yu  , b , and  X . For  the  more 
highly  excited  systems,  however,  the  influence  of  y , 6 and  X on 
Q.  is  appreciable;  consequently,  these  terms  must  be  included. 

The  partition  function  is  proportional  to  (Cf.  Equation 

(44) );  therefore,  it  is  of  interest  to  consider  the  evaluation  of  the 

|1Q\ 

infinite  series  (45)i  The  Euler- Mac laur in  formula ' ' may  be  ap- 

plied for  this  purpose.  If  (45)  is  written  in  the  form  of  Equation  (42), 
and  the  infinite  series  in  the  expression  approximated  oy  the  first  two 
terms  of  the  formula,  then 

J\ 


S “ Sli  + 2 

r 4 a 


(50) 
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The  function 


cj(n)  is 


tx/h2-/S(nSA,n4+^n2  + r) 

C|(n)^20n  e =200  h(n) 


(51) 


The  problem  is  to  evaluate  the  integral  of  Equation  (50).  Unfortunately 
the  integral  of  <^(n)  even  in  its  present  form  is  an  unusually  complicated 
function  and  is  not  amenable  to  any  elementary  methods  of  analysis. 
Numerical  calculations  indicate,  however,  that  the  function  h(n)  de- 
fined by  Equation  (51)  may  be  approximated  in  the  neighborhood  of  some 
value  n-a^  by  a linear  function  (Cf.  Figure  2).  Therefore,  we  expand 
h(n)  in  a Taylor  series  about  the  point  n=a^  and  take  the  first  two 
terms. 

h(tj)=»  h(ai)  + h'(at)(y-at)  + ••  • ;52) 


where 


- /n2 


h(q,)=  e 


h(at) 


If  these  expressions  are  substituted  into  the  above  mentioned  integral, 
then 

00  nQvb 

29c/yThfa*)+ 

a a J 

The  upper  limit  o o has  been  replaced  by  which  denotes  the  zero 
of  the  approximation  (52).  The  integral  is  terminated  at  in  order 

to  avoid  the  negative  values  predicted  by  the  approximation. 


-32- 


lb 

v 


h(d«) 

h'(cg 


(53) 


The  substitution  of  this  relation  Into  the  integral  yields 


Jm  “ 23 sH*/]] 


hn  \.~A  f 1 Ti 


This  expression  can  be  simplified  by  a suitable  choice  of  which 
denotes  some  convenient  point  in  the  range  of  h(y)  wherein  the  approx* 
imation  (52)  is  valid.  In  this  region  the  function  0j  , which  may  be 
represented  by  , is  of  order  one.  Thus  we  are  interested  in 
values  of  n—  0(jb  5)  . Let  us  select,  therefore,  the  following  value 
for  a.  : ' 


°*3  (13p) 


i 

<6 


(54) 


This  assumption  will  be  applied  to  all  calculations  so  long  as  is 
sufficiently  small  to  yield  values  of  > a . If  the  value  computed 
from  Equation  (54)  is  less  than  o , then  the  value  a will  be  assumed 
instead. 

The  particular  merit  of  assumption  ( 54)  lies  with  the  fact 


B,  and  (13^8)  6 is  a good  approximation  to  the  maximum  point  of  this 
function.  If  Equation  (54)  yields  a value  for  a^>  G , then  is  in- 
deed a maximum  point  of  the  function  <^(n) , and  it  follows  by  differ- 
entiation of  Equation  (51)  that 

The  substitution  of  this  relation  into  Equation  (53)  yields 


•3j* 


a«~la. 


and 


'W 


hj  ^ ea^hfa^) 


a 


where 


e = e 


The  internal  partition  function  may  be  approximated  by  means 
of  the  above  relations.  Thus 


*1  = ^ 


W-AT  V"_  .1 3 


I 


(55) 


The  second  and  third  terms  represent  the  infinite  series  as  previously 
indicated. 

Figure  3 exhibits  the  accuracy  with  which  the  function  <^(n)  can 
be  approximated  for  both  small  and  large  values  of  according  to 
the  relations  developed  above.  For  small  jb  » the  approximation  pro* 
vide s an  adequate  match  to  the  true  function  in  the  neighborhood  of  the 
maximum  contribution  of  <^(n)  to  . For  larger  jb>  « the  comparison 

is  not  as  good.  In  fact,  because  of  the  selection  of  according  to 
Equation  (54)  and  the  assumption  of  zero  slope  for  <j(n)  at  Qt  , the  ap- 
proximate  function  has  neither  the  correct  slope  or  curvature.  This 
problem  was  investigated  in  more  detail  in  order  to  obtain  a better 
match.  An  improvement  upon  the  above  mentioned  method  was  at- 
tempted by  inserting  the  correct  value  of  the  slope  for  Cj(n)  at  n=  0^, 
and  by  considering  various  other  positions  for  . In  none  of  the 
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* 


caves  investigated,  however,  did  the  corrected  relations  yield  as 
good  an  estimate  of  the  integral  (C£.  Equation  (50) ),  as  the  pro- 
posed approximations.  It  appears  that  this  method  is  the  best  if 
only  the  linear  term  is  to  be  retained  in  the  expansion  of  h(y) . The 
accuracy  can  be  improved  however  if  the  curvature  term  is  included. 
For  the  present  analysis  the  two  term  expansion  will  b«  sufficiently 
accurate. 

The  actual  function  c^(n)  la  compared  graphically  in  Figure 
3 with  the  approximate  expression  involving  the  two-term  expansion 
of  h(n)  as  indicated  in  Equation  (52).  The  actual  sura  is  also  com- 
pared with  the  two-term  Euler-Maclaurin  expansion,  and  for  both 
values  of  temperature,  the  agreement  was  found  to  be  good.  Com- 
putations indicate  that  for  T=5,ooo  the  error  in  Q is  leas  than 

mt 

1 percent,  whereas  for  J=2.opoo  , the  error  is  less  than  9 percent 

The  approximation  to  given  by  Equation  (55)  may  be  ap- 
plied to  the  solution  of  the  simultaneous  Equations  (47).  The  cor- 
responding relations  are 


.■f +lVa)  + ea*h(°*)} 

HI 

iVtV' + €a»hw] 

ii 

o-OI  — 

A.  

^<%lV)+6a*9h(a*)| 

where  the  notation 

has  been  introduced.  The  solution  of  Equations  (56)  by  the  suggested 
iteration  procedure  is  straightforward  but  lengthy.  For  this  reason 
the  derivatives  indicated  in  Equation  (49)  have  been  computed  and 
summarized  in  Appendix  C. 

It  is  a relatively  minor  problem,  however,  to  investigate  the 
solution  of  these  equations  for  the  extreme  values  of  temperature. 

At  very  large  temperatures  the  major  contribution  to  the  internal 
partition  function  is  from  the  higher  energy  states.  Therefore,  the 
£ Q*  h (cf)  term  in  Equation  (55)  is  predominant,  and  a good  approxi- 
mation to  is. 


^im  Q = £G^  h(d#) 

T-*-oo 

Therefore,  by  neglecting  the  first  two  terms  in  each  of  the  Equations 
(56);  we  obtain  immediately  the  following  asymptotic  values  for  the 
three  moments; 


,u-3a*  S-5a,4  *=  a*  <58> 

We  are  concerned  here  with  the  case  for  . Thus  the  function 

<^(n)  is  of  the  type  indicated  in  Diagram  B,  and  a*  represents  in 
fact  the  maximum  point  of  the  function.  It  is  possible  to  derive  the 
appropriate  value  of  from  the  asymptotic  expressions  developed 
above.  The  condition  of  zero  slope  for  <^(n)  at  n=  <3^  , along  with 
the  values  given  in  Equations  (58),  yield  the  following  equation  in  a*  : 


O 6 2 

13/3 Qt  - c x/3at-  ot  -f-  * o 


(59) 


For  the  cases  in  question  a representative  set  of  the  reiati'/c  magni- 
tudes involved  are  o^a^O^O)  and  $ — 0(i0'4’).  Therefore,  one  can 


obtain  an  approximate  solution  to  Equation  (5?)  by  neglecting  the 
a*  and  terms,  The  corresponding  solution  is  the  relation 
previously  assumed  for  ot  and  presented  as  Equation  (54). 

The  investigation  of  the  solutions  for  the  simultaneous  Equa- 
tions (47)  in  the  limit  for  very  low  temperatures  also  yields  a set 
of  simple  expressions  for  the  moments.  For  this  case  the  major 
contribution  to  Q is  from  the  first  term  (ground  state).  There- 
fore,  if  we  neglect  all  others , then 


If  these  expressions  are  substituted  into  Equations  (47),  one  ob- 
tains readily  the  low  temperature  solutions  for  the  three  moments. 


P = 3vt  S-3^2  r=V3 


(60) 


These  two  sets  of  solutions.  Equations  (58)  and  (60)  can  be 
used  effectively  as  the  initial  estimates  jj0  , Sa  and  X0  which  are 
required  in  the  iteration  procedure  outlined  above.  F igurc  4 gives 
the  exact  solution  to  the  simultaneous  Equations  (56)  as  a function 
of  temperature  for  a given  pressure.  It  is  interesting  to  compare 
these  solutions  with  the  high  temperature  approximations  (broken 
line  curves)  which  have  also  been  included  and  to  note  the  asymp- 
totic behavior  of  the  two  sets  of  curves. 
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V.  ATOMIC  VOLUMES 
A.  Hydrogen-like  Model 

In  the  present  section  we  examine  a method  for  computing  the 
atomic  volumes  which  determine  the  weighting  function  according 
tc  Equation  (12).  For  the  purpose  of  computing  thermodynamic  func- 
tions. it  will  be  desirable  to  sacrifice  accuracy  for  the  sake  of  sim- 
plicity and  to  consider  approximate  methods.  In  order  to  obtain 
some  information  regarding  probable  errors  incurred  by  the  approxi- 
mations, the  familiar  method  of  screening  constants^’  ^ was  applied 
as  the  approximate  technique,  and  the  results  in  a few  particular  cases 
compared  with  more  accurate  solutions  obtained  by  the  Hartree  meth- 

/ A\ 

od.  'w'  The  present  discussion  is  confined  to  a review  of  the  approxi- 
mate method. 

The  application  of  the  method  of  screening  constants  to  the 
calculation  of  atomic  volumes  consists  of  treating  a given  atom  as 
an  hydrogen-like  atomic  structure.  This  approximation  would  ap- 
pear to  be  reasonable,  especially  for  the  higher  electronic  states, 
on  the  basis  of  the  argument  already  given  in  Section  III  B.  For  the 
lower  energy  states,  the  shielding  effect  of  the  outer  electrons  is 
not  as  effective;  hence,  it  is  necessary  to  apply  corrections  in  the 

computation  of  the  effective  core  charge  by  means  of  empirical 

14) 

rules.  Slater  suggests  a set  of  rules'  based  on  data  obtained  from 
ions,  x-ray  energy  levels,  etc.  The  total  shielding  effect  of  the 
outer  electrons  on  the  field  of  the  nucleus  is  determined  by  com- 
puting the  contributions  ot  each  electron  by  means  of  the  rules. 


The  method  is  readily  applied  to  volume  calculations , and  is  there- 
fore especially  suited  to  the  present  problem. 

It  was  previously  noted  that  the  term  values  for  the  hydrogen 
atom  can  be  approximated  by  the  Bohr  model. 


/\J 


Z2R0hc 

nz 


(25) 


where  Z , (the  atomic  number)  is  of  course  one.  If  the  many-electron 
atom  is  treated  as  an  hydrogen-like  structure,  it  is  necessary  to  in- 
troduce an  effective  atomic  number  Z ..  which  is  defined 

eff 

z„,  a(z-5)  if") 


where  S is  the  screening  constant  computed  by  the  method  indicated 
above.  The  analogous  energy  relation  (in  Rydberg  units)  to  expres- 
sion (25)  is,  therefore 


W. 


n t 


R.hc 


(Z-  5) 


ni 


(62) 


It  is  assumed  that  this  relation  will  be  valid  for  any  arbitrary  atom. 
Note  tliat  spectroscopic  data  yields  £nt  which,  along  with  the  com- 
puted values  of  5 , define  n#  . 

The  calculation  of  atomic  volumes  is  made  on  the  basis  of 
an  assumed  spherically  symmetric  configuration.  Thus  the  determi- 
nation of  the  average  value  of  the  atomic  radius  cubed  (P*)  is  re- 
quired. This  quantity  may  be  computed  readily  from  the  relation 

00  JI  2J\ 
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where  r , and  (p  denote  the  polar  coordinates,  and  ^ i3  the 
wave  function  of  the  system  in  question. 

In  the  case  of  the  hydrogen  atom  it  is  known  that  the  wave 

(19) 

function  is  given  by 


Tj,  =4(r,i3-,<f>)  - Rr<(D 


(63) 


S (l^cp)  is  a normallized  surface  harmonic  of  order  |m|  and 

degree  t , and  Rnl(p)  is  defined  a3 


R„((r>  = 


\naj 


3 3 

(n-i-i) 

2n[(n+t)l] 


h 2 


-r\/z  l 
e r]_ 


2H1 

l (n) 

n+l  <- 


L (n)  is  the  associated  Laguerre  polynomial  of  degree  (£-n) 

n+£  l 

and  order  Zi  1,  and  rj=2Z^/na0  with  00  « radius  of  the  first 
Bohr  orbit.  It  follows,  therefore,  that 


(r_ \3  = /n  f (n-j-e)l 
[?zb_n[(n  H)!] 


(64) 


The  integral  may  be  readily  evaluated  to  yield  (Cf.  Appendix  B) 


(rH-l  + 4-fc)! 


tr-i4) 


It  should  be  noted  that  this  sum  terminates  at  if  n-£  <5  . 

If  the  summation  is  expanded  out  and  the  resulting  expression  substi- 
tuted into  Equation  (64),  then 

°4„  - tetrj  1 " ' (l  n '1*  + -n4tt4 

"N 

I 36 (n-M)(i?-t-2)  16(n-l-i)-(n-£-3)  (n-M)- (n-J-4)  [ 

(n+£+4)(n+i+3)  (n+t+4)-"(n+I+2)+  (n+U4)- (ntt+i)  f 
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which  may  be  reduced  to 


fuf-  ''ffiJ  3/J.f  , 6/lf_  6/iV5+l  U/35+2S\ 

\Z  k[  ln)+H\nl  1 77/ V + *V  n(n/{  +n4/ \ + n2/ 

t 

Or,  if  the  quantity  in  the  brackets  is  denoted  by  , then 


V (65) 


(66) 


An  analogous  relation  may  be  assumed  for  the  nvmy -electron  atom 
by  replacing  the  (nzf Z ) factor  by  Equation  (62), 


(67) 


The  spherical  volume  b,  corresponding  to  this  value  of  is 

,68> 

where  the  more  general  notation  i for  the  quantum  energy  state  has 
been  introduced.  Xt  is  evident  from  Equation  (26)  that  in  computing 
the  internal  partition  function,  it  is  necessary  to  sum  the  product  of 

-a 

the  degeneracy  and  the  weighting  function  q e over  the  allowed 

'-'In 

values  of  l , The  computation  if  complicated  by  the  presence  of  the 

function  ){.  , which  is  defined  by  Equation  (66).  Figure  5 gives 

\ n 

some  indication  of  the  nature  of  this  function.  Note  that  1 < V.  < 35 

qln  & 

and  also  that 


- §(35+  %) 
W = 1 + n + 4? 
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A general  approximation  can  be  obtained  from  Equation  (65)  for  small 


In  spite  of  these  simplifications  the  summation  over  l proved  to  be 
a formidable  problem  and  no  simple  solution  was  found.  In  view  of 
thi3  fact,  it  seemed  advisable  to  si  oiplify  the  entire  problem  by  as- 
sijnir.g  to  J(^n  some  average  value  )(0  . Equation  (34)  is  the  cor- 
responding simplification  of  the  general  expression  (68). 

The  selection  of  )l0  for  a given  system  car.  be  accomplished 
by  an  iteration  process.  Starting  with  an  initial  estimate,  one  can 
compute  by  first  solving  for  the  moments  JJ  , £ and  X . An 

average  value  of  V,  can  then  be  defined  for  each  value  of  n . If 

m 

the  computed  average  compares  favorably  with  the  assumed  value 
in  the  region  of  n wherein  cj^(n)  makes  the  largest  contribution  to 

q , then  the  initial  estimate  was  rwa  adequate  one.  If  not,  then  the 

6 

process  can  be  repeated.  A suitable  average  )j.  for  given  n can 


be  defined  as 


-S^V/inV  + irfys) 


n 1-0 


where 


Ss4 
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this  function  as  the  assume.!  value  in  varied  from  K>  * 2 to  V0*  2.  6. 
The  influence  oi  on  ^ is  also  a relatively  small  effect.  Sample 


calculations  indicate  that  for  « 0(10  ),  a 25u/o  increase  in 


causes  about  a 5 /o  decrease  in 


% 


B.  Calculation  of  Atomic  Volumes 

Atomic  volumes  for  any  quantum  energy  state  i-  may  be  com- 
puted for  a spherical  model  (baaed  on  the  assumption  of  a hydrogen- 
like structure)  by  using  Equation  (67).  The  quantity  V.  is  computed 
from  the  known  spectroscopic  data  of  the  atom  (assumed  hydrogen- 
like for  the  higher  states),  and  the  screening  constant  S , by  tha 
method  previously  mentioned.  It  is  of  interest  to  determine  the  value 
of  for  the  higher  energy  states  according  to  the  hydrogen- 

like approximations.  In  these  cases,  (Z~S)-*-l  and  d = — 2 from 

In  r, 

Equation  (62).  It  is  evident  from  Equation  (67)  that 


V.  = V = n 2 

In.  n 

and 

b.  b,  = ( 4JTg.»3 ) V n6 

1 \ 3 / fn 


i 69 ) 


The  calculation  of  these  volumes  may  also  be  carried  cut  by 
applying  the  more  accurate  method  of  Ilartree,  usually  referred  to 
as  the  method  of  self-consistent  fields.  The  procedure  is  essentially 
a numerical  one  and  relics  primarily  on  the  representation  of  a com- 
plex atomic  structure  by  the  superposition  of  single -electron  wave 
functions.  A numerical  representation  of  the  wave  function  for  the 
entire  system  is  obtained  by  trial  and  error  estimates  on  these 
single-electron  functions.  The  applicability  of  the  method  to  the 


-43- 


present  problem  was  Investigated  by  ..ctually  carrying  out  volume 
computations  (r  )3  for  the  potassium  atom,  rise  semi -empirical 
wave  functions  for  this  atom  were  supplied  by  D.  S.  Viilars  of  the 
U.  S.  Naval  Ordnance  Test  Station,  Inyokern.  These  data  were 
tabulated  for  the  radial  wave  function  P(r)  as  a function  of  the  rad- 

in 

iai  distance  from  the  center  of  the  nucleus.  The  function  P (r)  is 

In 

d‘  lined 

P (r)  = p R (r) 

In  |r> 

in  accordance  with  the  expression  for  Q given  by  Equation  (63). 

The  function  (P)3  was  readily  obtained  by  the  evaluation  of  the  fol- 
lowing integral  (Cf.  previous  definition  of  r 3 ) using  numerical  meth- 
ods of  integration: 


Oo 


o 


Sample  computations  of  (P)3  for  various  energy  states  are  presented 
in  Table  I and  compared  with  the  results  obtained  by  the  method  of 
screening  constants.  The  correlation  is  well  within  the  limits  of  ac- 
ceptable engineering  accuracy. 
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VI.  THERMODYNAMIC  FUNCTIONS 


A.  Fundamental  Relations 

The  following  well-known  relations*  define  the  contribution 
to  the  thermodynamic  functions  by  the  internal  energy  states.  In 
the  present  case  Q is  due  to  electronic  excitation. 

tat 


_EL  = _fnQ 
kX  int 


JL  = id  = _L_  f j ^9{oi) 
kJ  kT  Q V dl  J 

int 


(70) 


where  F is  the  Gibbs  free  energy 


F - H-  TS 


(71) 


The  heat  capacity  C is 


C 

k. 


(72) 


The  evaluation  of  these  functions  requires  explicit  relations  for  Q.  ^ 
and  its  first  two  temperature  derivatives.  These  derivatives  may  be 
computed  in  general  form  directly  from  the  Fermi  Equation. 


B.  Thermodynamic  Functions 

Good  approximations  for  the  thermodynamic  functions  can  be 
obtained  by  using  Equation  (55)  for  and  performing  the  neces- 

sary differentiations  directly  on  this  expression.  In  carrying  out 
thin  scheme,  it  is  convenient  to  separate  Q into  two  parts  cor- 

int 

responding  to  the  notation  of  Equation  (26).  Thus  define 

* Cf.  Reference  10,  p,  442. 
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0 .=  Q(°\  Q(0 

int  ^ 

wh <9 re  according  to  Equation  (55) 


q1.wwiAt£ 

r 1 J 


Q - Q0  + Q*  ~ IS(o)e  " + €o;h(at)e 


-vviAT  3,  / , ^/o 

fhe  first  derivative  of  Q .n{,  with  respect  to  T is 


T £9 h*  - T^'-f  T&Qo'  + T £Q* 

aT  a~r  ^ aT 


.<»> 


(U 


(73) 


where 


T == 

' aT 


^ - 
CO 

T 3S#- 


3T 


(>) 

Q0  A(a) 

q“a  (0*) 


and 


I 


a(q)  = S'  + fj  {a^  + ^ + ha  +?) 

kJ 


Wr  ~ W,- 
a 1 a2 


(1-jVt  )2 


The  second  derivative  of  Q)-n£  is 


3t2  aT2  + aT2  at? 


1741 
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where 


T'&r.  e~w,AT. 

3T2 


I f 

J 


2 \zs-\i)  W 


A (a)  -m-  IA  (aSu04^o!tri 


1 

j 


t‘£&  - Q^A^-?r^-(T^_a(a>o;+H£^)  ■ 


The  substitution  of  these  derivatives  into  the  expressions 
for  the  thermodynamic  functions.  Equations  (70)  to  (72),  is  a straight* 
forward  procedure. 

C.  Equation  of  State 

An  equation  of  state  for  a gaseous  system  of  neutral  atoms 
with  excited  electronic  states  can  be  derived  in  terms  of  the  results 
which  have  already  been  obtained  in  the  previous  sections.  In  par- 
ticular, it  can  be  shown  that  the  covoiurne  correction  b given  in 
Equation  (1)  and  defined  in  Equation  (4)  may  be  expressed  as  a func- 
tioii  of  the  moments  of  the  internal  partition  function  (Cf„  Equations 

(38))  and  the  state  parameter  /i  . 

I 

Consider,  therefore,  the  general  relation  (1).  This  expres- 
sion may  be  simplified  considerably  by  introducing  the  relations 
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(27),  (33)  to  (35)  and  (44),  Thus 


where 


b ,= 


2JTaoK, 


(75) 


I 


If  the  above  result  is  expanded 


s “ + 3V‘  phrfftl 


Obviously 


0 or- 

jp_  = i yv-v:3+  ^2+  -l  r; 

b„N  cjfclH  L+/  1 L+  j 


And  so 


JL  « 2(t+  4m) 

bN  v s' 


(76) 


The  corresponding  equation  of  state  is  readily  obtained  by 
substituting  this  relation  into  Equation  (1),  which,  in  standard  form 
is 


-PV  . 1+  bP  ^ , b 
NJcT  NkJ  r V 


(77) 


If  the  definition  of  /b  is  applied  along  with  Equation  (76) : the  above 
equation  reduces  to 

PX  « l + zp(7t-f) 


(78) 


It  is  of  interest  to  examine  the  behavior  of  this  function  in 
tiie  limit  for  large  temperature.  According  to  the  analysis  of  Part  IV 
B above,  the  high  temperature  limits  to  the  three  momenta  are 
given  by  Equations  (58).  These  results,  along  with  the  approxima- 
tion (54),  yield 


lim  (Py)  = 

J-+.00  ' RT/ 


u 


8 

13 


(79) 


Some  sample  computations  of  the  "compressibility  factor"  given  in 
(78)  are  presented  in  Figures  7,  8 and  9. 

The  equation  of  state  (78)  may  be  corrected  to  account  for 
the  presence  of  ions  in  the  gas  mixture  which  appear  as  products  of 
the  reaction  (28).  In  making  the  actual  calculation,  it  will  be  more 
convenient  to  use  a slightly  different  definition  than  indicated  in  liqua- 
tion (77).  Consider  therefore  the  expression 


PV  _ j , _Pb_ 
7\lcT  T\kT 


(SO) 


where  ^ denotes  the  total  number  of  particles  at  equilibrium  accord- 
ing to  the  definition  of  Part  III-C.  The  substitution  of  Equations  (29) 


into  (30)  yields 


PV „ 1+  Pb 

Nkj”  (J  + ) T\’<T 

which  is,  in  the  notation  of  Equation  (73) 


PV 

NkJ 


(ai) 


Note  that  the  limit  of  the  ionization  factor  for  large  T is  2. 


-49- 


D.  Discussion 

The  present  analysis  has  achieved  three  objectives.  First, 
a convergent  series  expression  for  the  internal  partition  function  of 
un-ionized  atoms  has  been  derived.  Second,  this  result  provides  a 
method  for  computing  the  electronic  contribution  to  the  thermodynamic 
functions.  And  finally,  the  range  of  applicability  of  both  the  customary 
method  and  of  the  covolume  treatment  has  been  determined.  The  re- 
sults of  the  analysis  show  decisively  that  although  the  covolume  cor- 
rections are  the  major  deviations  from  ideal  gas  behavior  at  the  lower 
temperatures,  these  effects  may  be  neglected  in  making  engineering 
calculations.  Furthermore  it  is  found  that  at  the  higher  temperatures, 
when  these  deviations  become  significant,  it  is  necessary  to  include 
also  the  effects  of  ionization.  Indeed,  actual  computations  performed 
for  lithium  and  nitrogen  (Cf.  Figures  7,  8 and  9)  indicate  that  the 
non-ideal  behavior  of  the  gas  at  the  higher  temperatures  is  due  large- 
ly to  this  effect,  and  that  the  covolume  correction  to  the  equation  of 
state  due  to  the  expanding  volumes  of  the  excited  atoms  is  in  com- 
parison a smaller  deviation.  Therefore,  if  these  results  are  at  all 
representative  of  the  general  behavior  of  all  atoms,  then  it  is  appar- 
ent that  a realistic  derivation  of  the  thermodynamic  functions  of  gases 
at  high  temperatures  must  necessarily  take  into  account  the  presence 
of  ions  and  free  electrons  in  the  mixture. 

The  extension  of  the  present  treatment  to  include  the  effect 
of  ions  and  free  electrons  in  the  mixture  is  a formidable  problem. 

The  most  serious  complication  which  arises  in  tills  calculation  is  the 
specification  of  the  interaction  between  the  various  types  of  particles. 
Although  the  method  of  excluded  volumes  has  been  applied  to  approx  “ 
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im&te  the  interaction  of  neutral  atoms  in  excited  electronic  states, 
the  application  of  this  method  is  not  generally  recommended  for  sys- 
tems which  contain  charged  particles.  An  honest  calculation  of  the 
thermodynamic  functions  of  a mixture  of  atoms,  ions  and  electrons 
must  be  based  on  a more  realistic  description  of  the  Interaction  po- 
tentials; but,  these  potentials  are  generally  not  known. 

When  ionization  is  of  the  order  of  10  percent  (at  temperatures 
of  about  10,  000°K)  charged  particles  will  be  Involved  in  relatively  few 
encounters,  and  it  is  expected  that  the  interaction  terms  corresponding 
to  binary  collisions  between  charged  particles  and  neutral  atoms  will 
yield  small  corrections  on  the  thermodynamic  functions.  It  is  suggested 
therefore  that  as  a first  approximation  these  interaction  effects  may  be 
neglected  entirely,  and  that  the  mixture  may  be  treated  by  the  familiar 
methods  of  classical  solution  thermodynamics,  In  these  calculations 
one  may  compute  separately  the  contributions  to  the  thermodynamic 
functions  of  the  atoms , and  of  the  ions  and  free  electrons  and  then  apply 
the  necessary  mixing  terms  to  determine  the  properties  of  the  complete 
system.  The  charged  particles  i n the  mixture  may  be  treated  by  the 
Debye-Huckel  theory and  the  methods  outlined  in  the  present  study 
are  directly  applicable  to  the  neutral  atoms . These  approximations  are 
admittedly  crude,  but  an  attempt  to  achieve  greater  accuracy  would 
require  much  more  realistic  descriptions  of  the  interactions  for  the 
various  atomic  particles. 

The  relations  derived  in  the  present  analysis  have  been  applied 
to  the  computation  of  the  thermodynamic  functions  of  lithium  1 
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for  various  temperatures  and  pressure.  In  particular,  the  contri- 
butions from  the  electronic  states  to  the  internal  energy,  the  entropy 
and  the  heat  capacity  have  been  plotted  as  functions  of  temperature 
in  Figures  10  to  12.  It  should  be  noted  that  these  functions  represent 
the  contributions  of  only  the  neutral  atoms.  Broken-line  curves 
have  been  drawn  for  the  higher  temperature  range  of  the  functions 
because  of  the  uncertainty  in  the  method  of  computing  atomic  vol- 
umes for  the  more  highly  excited  atoms.  TheBe  volumes  are  rela- 
tively small  for  the  lower  energy  states,  and  large  errors  can  be 
tolerated  in  the  calculation.  But,  the  influence  of  the  atomic  vol- 
umes (i.  c.  interactions)  on  the  thermodynamic  functions  become 
more  pronounced  at  the  higher  temperatures:  thus,  more  realistic 
interaction  potentials  should  be  introduced  for  atoms  in  the  highly 
excited  electronic  states. 
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TABLE  I.  ATOMIC  VOLUME  CALCULATIONS 
FOR  POTASSIUM  L 

E“"  ' /z-jA  T y [75  I" (tf 

emission  \JT  ) i("  (Sister's  \^"j 

electron  method)  (Hcrtree) 

3d  541  2.33  1260  102? 

5-4  516  4.  50  2320  2366 

5d  11,750  3.61  42.400  40,326 

7-4-  12,500  4.44  55,400  55,900 

114-  479,000  4.37  2.090,000  2.062,665 


TABLE  H.  THERMODYNAMIC  FUNCTIONS  OF  UTHIUM  L 


a 
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Temperoture  T (0.K,) 
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<0 
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2.3  • 10 
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Temperature  T (°KL) 
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7- 
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erature  T (X) 
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APPENDIX  A 

Classical  Solution  Thermodynamics 
1 . The  Chemical  Potential 

For  real  gases,  the  fugacity  fL  of  component  i is  proportional 
to  the  activity  \ . 

■ const,  (A»l) 

The  activity  , however,  is  related  to  the  chemical  potential  ^ by 

e/J‘/llT  (A-2) 

If  we  substitute  Equation  (A-2)  in  (A-l)  and  solve  for  ^U;  , we  obtain 

w=tTH^) 

Let  the  properties  of  this  component  in  the  pure  state  ( ) be  used  as 
a reference,  then 


^hSt) 


and  it  follows  that 


rt-K  + iTHv* () 

But  the  mole  fraction  is  given  by 


Ni 

N 


l*  4K  mm  mm  rm  ♦ 

UQAAWC* 


fi 

ft" 


(A-3) 


2.  The  Entropy 

The  differential  equation  defining  Gibb's  free  energy  F is 

ciF  = -SdT  + VdP  4.  £ ft  city 

c 


(A -4) 
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but 


F“  H/'.N i 

i. 


(A-5) 


00  that 


Jr 

VJ  I 


N ' I . ,-JK  I.  \ l M I . .. 

= n + ZV\un 

t i 

If  we  equate  this  expression  to  (A-4),  we  obtain  the  Gibbs -Duhem 


relation 

_ SdT-u  VdP-  « O 

t 

For  a constant  pressure  process  dP=o  and 
SdT+  ° 


(A-6) 


(A-7) 


The  entropy  may  also  be  defined  by  Euler's  relation  for  extensive 
properties 


PTM; 


(A-8) 


If  Equation  (A-8)  is  substituted  in  (A-7) 
But  is  arbitrary,  therefore 

/ ^ ! I \ 


iar/  UnJ 

PNj  l 


PTNj 


(A-9) 


We  can  obtain  the  temperature  derivative  of  jj i from  Equation  (A-3) 
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L In 

N 


But  for  a pure  component  Equation  (A-7)  yields 


(A-10) 


Si  _ _ (A»li) 

Ni  '5T  'pNj 

And  from  (A-9),  using  (A-10)  and  (A-ll) 

_ L In  Hi- 

V3N^pTN.  Nj,  N 

the  substitution  of  this  expression  in  Equation  (A- 8)  gives 

S = 2]($L“  (A-12) 

l 

% 

3.  The  Volume 

For  a constant  temperature  process,  the  Gibbs-Duhem  rela- 
tion (A- 6)  yields 


V- 


’ K I / 


wp 


T 


The  Euler  relation  for  the  volume  is 


(A=13) 


V 


- 


UsnVptk 


'i 


(A-14) 


For  a pure  components  Equation  (A-13)  reduces  to 
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Vi=lhd\  (A  *15) 

Nt  \dP/T 


By  using  the  pressure  derivative  of  ^ from  Equation  (A*3)  we  obtain 


(A-16) 


where  Equation  {A-15)  has  been  introduced.  The  substitution  of  this 
egression  into  Equation  (A-13)  yields 


v-Ev- 

i 

and 


(A-17) 


(A-18) 


4.  The  Internal  Energy 

The  Gibbs  free  energy  F may  also  be  defined  by 


F as  H - TS  - E +-  PV  ~ T 3 {A *19) 

If  we  introduce  Equation  (A*5)  into  this  expression  and  solve  for  E , 
the  internal  energy. 

E » T5  - PV  + E Nj  u 

i ' 

The  application  of  the  relations  (A *12)  and  (A *18)  yields 

E “ Z,'  |T(S[  - k-N;  - PV,"  +/liNtj 


(A-20) 


The  expression  for  from  (A-3)  may  now  be  introduced 

E“Z}TSt'-  pV+  nlK 


where  the  quantity  in  the  bracket  is  evidently  the  internal  energy  of 

EO 

I , hence 
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APPENDIX  B 

1.  The  Associated  I^aguerre  Polynomial 


L"  “ d^Ln« 


(B-l) 


where 


] =|o^p 


(B-2) 


If  Equation  (B-2)  is  differentiated  k -times 

u?s*p\zm. l . [-z_ik  = y 

i-t  1 Li-  Z.J  L'-2U  fyk  A! 


(B-3) 


Define 


CX* 

V = /e  -X  LnWL'*'” 


where:  k n.ri.rn  are  integers 

* r~  * • w 


(B-4) 


Multiply  Equation  (13-3)  by  another  L,  a guerre  function  evaluated  at  Z=  Z2 

2 4nfr  ^ = (z'  z'>k  [('  -x-rw  [-  - ffe] 


Multiply  each  side  by  q ^ ^ P J 


OO 

2 


a ,« 


and  integrate 


DO 

-rK~L  / k-rO-i 


={Z'^]  J*  ei?tx~ &r  &Jdx 


1 
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For  ? ^ o , this  expression  readily  reduces  to 


OO 


V1 

/-J 


1 • 
*>/' 


l 


~k 


* M 

^ !/jl 


(AX.)k[0-^,Xi-z;)]P"‘(l<+P-l)l, 
(J  -?tZ*)k+p 


(B-5) 


The  denominator  of  Equation  (B-5)  may  be  expanded  in  a Taylor  series 
about  the  point  %\%z~  ° 


(1-H,zt)'k'p  _ V (ktp+a-  0!  (z, Za)'*  (b-6) 

^=o  (k+p-i)U| 

The  term  Zz)J  ^ J may  also  be  expanded;  thus  the  expansion 

of  each  factor  about  the  point  z = o yields 


[0-z,)(i-z!)]P  1 


. y (p-0!(-z,)^ 


f<  (p-QK-z.)* 

t=>  t!(p-i-l)i 


(B-7) 


The  substitution  of  Equations  (B-6,  7)  into  Equation  (B-5)  gives  * 


V~k 


Z’i 


M 


Al 


r[- 


rO 


= A 


dt—  o 


OO 

/ 


5 H) 

d)  tl  f p-i-d)  j (p-i-t)  i (A  k ) I 


# id  .Ml 


where  V * A 

An  expression  for  the  function  I is  obtained  by  setting 
jj  = Am  and  equating  the  coefficients  of  the  like  powers  of  (z,z2) 
from  the  two  sides  of  the  abuve  eqr  uion.  If  it  is  noted  that  only  like 
values  of  d and  t yield  terms  in  (z(Z2)m,  then  Equation  (B-8)  re- 
duces to 
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y1  (Z,Z.)  1m, i 

itl  (ml)* 


f y [»-)l]^».)tw(P*v-»)! 

<l,t  = O -fc" (S  [t ! ( p-t-i  )ij2  (V-  Is  ) ! 


Y1  [(p-Oilfo+^-Q! ' fz, zz)v 
v-k  (V-U!  U(p~l>!]2 


, (z.z2)^2  ( , foz,)Vtp~'~[ 

[fp-2)!jz  T[(p-3)l£!]2  pP-OI]  ‘ J 


Note  that  the  summation  over  x terminates  when  t = p-l  • By  equating 
the  coefficients  of  (Z,Z2)  and  using  1^4-1 J=m, 


l 


nr^m 


P ~~i  * 

[ml(p-l)  \]ZJ]  (p  + m-l-^)l 

i>=o  [l>j(p-i-'i?)!]2(m-lc-^)l 


(B-9) 


The  function  1^  m for  various  values  of  p can  be  obtained  readily. 


m.m 


M? 

(m-k)| 


(B-lUJ 


D=  2 

I 

= (mlf  Y _ (m Y?(Zm-k+l) 

V-o  [(l-TJ)jtf  |]z(m-'iJ-k)|  {.'n-kjj 
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p«  3 


T 

•‘■rn)rr( 


2 


(m-T>+2)  | 

[(2-'P)!fl*]2(m-tf-k)l 


0=5 


* -MI0  /fern2-  6mk  + k%  <onri-3k+  2 ) 
(m-k)P  ' 

W,  = ^U!f  y1  _ fcL±tll[ 

v-o  [w-wif  (m-i/-j<)! 


(5-12) 


(B-13) 


In  the  present  analysis  it  will  be  of  interest  to  define 

m •=  r.  i l 

k.  = <L  C + i 

where  n and  (.  are  the  principal  and  azimuthal  quantum  numbers  re- 
spectively. Equation  (B-13)  yields 


w,--!><>!4!rx:  (n+j~aP+4)! 


(B-14) 


«► 

k> 


At. 
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APPENDiX  C 


The  present  section  summarizes  the  important  relations  re- 
quired for  the  solution  of  the  simultaneous  equations  (48}  using  the 
approximations  (49).  The  derivatives  of  computed  from  the  rela- 
tion (55)  are 


+Hw+€q?h(<v] 

U—  nft[f +lV,+eq*5hM 
$ r (4r)£{f+lVa)+£a*hM 

The  corresponding  derivatives  of  <p>  , % and  are 

ff-  S+  I$r  (3f Vf V I-V<3a>) + 

|f=  A (a 

i 

ip - +n<#+  v )J  -/[f  +^lVf + *n] 

+mm)(5a‘-M)+  {tghtyEpMxf-M)  j 

If  = 1>+  rpr  {3<f _ *f + f 1 + eaf 1 %>(*$  - s;J 

(T^-{3[lc',+^^wV?WJ+^^cs03-«s£fV^V<i^  <‘)J 

+ £>a4-S  ) + €^fh(cy  E(<y(iq^-5)j 
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where  the  function  f(0)  ie  defined  by 

Ra)=  i +^(aV^a4  + ^a2) 

In  carrying  out  the  actual  computations  for  the  three  moments 
by  means  of  the  above  relations  and  Equations  (48)  and  (49),  a slight 
simplification  can  be  accomplished  by  neglecting  the  weighting  func- 
tion in  the  lower  energy  states.  This  is  to  say,  one  can  assume  an 
approximation  for  the  sums  (Cf.  Equation  (57) ) which  is  inde- 

pendent of  B . Thus 


This  approximation  is  especially  good  for  the  lower  temperature  e»ys- 
tems  ( (6>10  ^).  However,  it  is  also  applicable  to  the  smaller  valuer. 
oi  [5  . It  is  applicable  in  the  first  case  because  the  entire  weighting 
function  £?  1 . Although  this  is  not  the  i.ue  for  ^ < 10  it  may 

still  be  vsed  because  the  contribution  to  of  the  terms  j = i to  \ is 
not  the  major  contribution.  \3  previously  discussed,  the  largest  con- 
tribution stems  from  the  higher  energy  states.  Thus,  relatively  large 
errors  can  be  tolerated  for  the  first  ^ terms  of  (1  . 


